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ABSTRACT 



Context. Hydrodynamical simulations of two giant planets embedded in a gaseous disk have shown that in case of a smooth convergent 

migration they end up trapped into a mean motion resonance. These findings have led to the conviction that the onset of dynamical 

instability causing close encounters between the planets can occur only after the dissipation of the gas when the eccentricity damping 

is over. 

Aims. We show that a system of three giant planets may undergo planet-planet scattering when the gaseous disk, with density values 

comparable to that of the Minimum Mass Solar Nebula, is still interacting with the planets. 

Methods. The hydrodynamical code FARG0-2D-1D is used to model the evolution of the disk and planets, modified to properly 

handle close encounters between the massive bodies. 

Results. Our simulations predict a variety of different outcomes of the scattering phase which includes orbital exchange, planet 

merging and scattering of a planet in a hyperbolic orbit. 

Conclusions. This implies that the final fate of a multiplanet system under the action of the disk torques is not necessarily a packed 

resonant configuration. 
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1. Introduction 

A multi-planet system made of Jupiter-sized bodies may reach 
completion when the planets are still under the influence of 
the disk. In this scenario the mutual gravitational perturbations 
compete with the gravitational torques exerted on the planets 
by the disk. These torques cause planet migration which is a 
strong source of mobility to the system and possibly causes 
the planets to come closer to each other. The mutual approach 
is supposed to end when a pair of planets is captured in a 
mean motion resonance (e.g. the 5:3, 3:2, 3:1 or 2:1). The most 
extensively studied scenario is that of two iovian mas s plan- 
ets embedded in a viscously e yolviii g disk. Kiev et al. ( 2004 
l2005h. iMasset and SneUgrov3 (l200lh and iPierens and NelsonI 
(|2008|) have shown via self-consistent disk-planet hydrodynam- 
ical simulations that the two planets become trapped into a res- 
onance and create a cavity encircling their orbits. This resonant 
configuration appears stable in these simulations which argues 
in favor of a scenario where the resonance locking should be 
broken only after gas dissipation possibly by the end of the gas 
damping or by planetesimal scattering. 

Damped N-body models have also been used to study the 
evolution of two planet systems where the disk- planet interac- 
tion is modeled as an external dissipative force, l ^e and Peala 
(2002), Beauge et al. (2006) and Lee et al. (2009) have shown 
that the resonance trapping is permanent if a strong damping of 
eccentricity is adopted. Even turbulence appears not to prevent 
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the resonance capture (iRein and Papaloizoul l2009h in this sce- 
nario. The escape from resonance can occur only if the eccentric- 
ity damping term used in the model is relatively low and, after 
the resonance breaking, planets may have close encounters. This 
seco nd approach chara c terize s the outcome of the simulations by 
Mo orhead and Adamsl (l2005h . where two planets may undergo 
planet-planet scattering after the escape from resonance trap- 
ping. However, hydrodynamical simulations have so far never 
shown such a behavior with the planets assumed to remain in 
resonance until the dissipation of the gas. The gas damping ap- 
pears to be dominating over the mutual gravitational perturba- 
tions even for Jupiter size planets, which supports the idea that 
the preferred outcome of the evolution of a multiplanet system 
under the tidal action of the disk is a packed resonant system. 

A system of three giant planets has proved to be dynam- 
icall y more prone to instabihty and to develop crossing or- 
bits (Weidens chilhng and Marz aril [1996: Chatte riee et alil2008h 
compared to the case of two planets. The mutual gravitational 
perturbations are stronger and may more easily overcome the 
disk-damping effect. The observed systems Ups And, 55 Cnc, 
HIP 1 48 1 0, HD74 1 56, HD 1 6069 1 and HD 37 1 24 also show that 
it is indeed possible to form more than two giant planets around 
solar type stars. For this reason we selected a scenario where 
three giant planets with masses close to that of Jupiter form in the 
disk and evolve under planet-disk and planet-planet interactions. 
We will show that the chaotic evolution, which is characterized 
by mutual close encounters between (giant) planets, is favored 
by planet migration. Our hydrodynamical simulations show that 
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most of the time the planets are quickly captured in a three-body 
resonance during their migration, despite their rapid drift rate, 
but this may become unstable. The growth in eccentricity leads 
the planets into crossing orbits, and they begin to evolve under 
the effects of mutual approaches. The following chaotic phase 
ends up either with the ejection of a planet, planet-merging or 
orbital reconfiguration. 



2. Numerical set-up 

The evolution of the planets and disk was modele d with the hy- 
drodynamical code FARGO-2D-1D JCrida et al.L |2007) where 
the numerical integrator computing the planet orbits was up- 
dated with a variable stepsize so that mutual close encounters 
between the planets could be properly modeled. The 2D grid, 
which resolves the local disk-planet interactions, is surrounded 
by a ID (non azimuthally resolv ed) grid that allow ed us to simu- 
late the disk's global evolution ( Crida et al.L |2007|) . The 2D grid 
used in the simulations to represent the disk has a logarithmic 
radial spacing extending from r=0.5 to r=15 AU. It is character- 
ized by A^r = 512 elementary rings, and A^s = 256 sectors. The 
ID grid spans a radial distance from 15 AU outwards till 50 AU. 
The disk-aspect ratio (H/r) was set uniformly to 5% and con- 
stant in time because the disk was assumed to be locally isother- 
mal. The gas kinematic viscosity is v = 10"^ in code-normalized 
units, which corresponds to an alpha parameter of 4 x 10 "*. The 
gas density profile declines as r"''^ and the value Eq at 1 AU 
was varied in the simulations, ranging from 0.25 to 1 Emmsn, the 
density at 1 AU predicted for the Minimum Mass Solar Nebula. 
In the simulations presented here the three planets have masses 
equal to 1.5, 0.9 and 1.2 Jupiter masses, while the initial semi- 
major axes are selected to be between 2 and 10 AU. They are all 
started on circular orbits and do not accrete mass from the disk. 



3. Results 

Figure 1 illustrates the different behaviors we observed in our 
simulations. After an initial fast migration, the three giant plan- 
ets become locked into a three-body mean motion resonance. 
However, the resonant trapping is not stable and the growth in 
eccentricity leads to planet-crossing orbits, which in turn start 
the chaotic phase. In the top plots we show the semi-major 
axis and eccentricity evolution of the three planets initially set 
on orbits with oi - 3.0AU,fl2 - 6.7AU,fl3 = 10.5AU with 
So - 0.4Smmsn- After repeated mutual encounters, an orbital ex- 
change occurs in which innermost planet becomes the outermost 
one. In the middle plots, a case of planet merging is displayed, 
where the two surviving planets are left on mildly eccentric or- 
bits. In this model the initial disk density is So = 0.65I;mmsn and 
the planets are started at ai = 2.0AU, fl2 - 6.5AU, ay, = lO.OAU. 
Very large eccentricities are achieved during the chaotic phase, 
but they are partly damped by the disk after the merging event. 
It is interesting to note that in this second case, after the planet- 
planet scattering, the two planets left on eccentric orbits are sur- 
rounded with eccentric gaps. This is illustrated in Fig. 2. The 
bottom plots refer to a simulation where two planets merge and 
are ejected out of the system on a hyperbolic orbit (Eo - ^mmsn, 
ai = 3.0AU,a2 = 5.5AU,a3 = lO.OAU). 

Different set ups were tested to check the robustness of our 
findings. In one series of simulations, the mass of each planet 
grows from to the final value on a timescale ranging from 10 
to 1000 yrs, the latter being comparable with the rapid gas-infall 
timescale typical of the final stage of a giant planet growth. In 
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Fig. 2. Disk surface density distribution after the planet-planet 
scattering phase shown in Fig. 1 middle plots. The two planets 
have developed eccentric gaps. 



another series, planets are first held on fixed circular orbits un- 
til they form gaps, after what they are let to evolve under disk- 
planet and planet-planet interactions. Very similar behaviors are 
observed in all these cases. The mutual gravitational interac- 
tions within a system of three planets force some eccentricity, 
which is not fully damped by the disk. In addition, the wakes 
excited by each planet tamper with the gap-cleaning process of 
the other planets. The two outermost planets in our simulations 
do not undergo slow type II migration. They experience strong 
torques from their co-orbital region w hich is not fully cleared . 
The fast migration (runaway migration, iMasset and Papaloizoul 
2003) triggered by corotation torques leads the planets close 
enough to be trapped into mean-motion resonance. But this is 
not stable and planet-planet scattering follows. 



4. Discussion and conclusions 

Although there are known two-planet systems in simple mean- 
motion resonances, there are also many systems that are either 
multiple and non-resonant, or that display evidence for scatter- 
ing that would be suppressed if resonant capture was a ubiq- 
uitous process. We show that resonant capture is more easily 
evaded in a three-planet as opposed to a two-planet system, 
leading to a phase of planet-planet scattering. A major impli- 
cation of planet-planet scattering in the gas disk is the signif- 
icant re-arrangement of planet orbits and eccentricity excita- 
tion, which may affect the subsequent evolution of the plan- 
ets. As a consequence, the final mass distribution of planets 
as a function of the semi-major axis cannot be representative 
of the formation process. In addition, a closely packed plane- 
tary system where the multiple planets are all locked in mean- 
motion resonances is not the preferred outcome of planet migra- 
tion, as previously argued. The chaotic phase may lead to dif- 
ferent orbital configurations where the planets can be well sep- 
arated into semi-major axis and stable even after the gas dis- 
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sipation. Lower eccentricities are expected as the final state of 
the planet evolution because damping can occur after the dy- 
namically active phase (Moeckel et al. , 2008). However, it may 
not be the case for corotati o n resonanc e saturatio n. Acc ording to 
iGoldreich and Saril (|2003|) : lOgilvie and Lubowl (1200 3') a planet 
may experience a net growth of eccentricity via disk-planet in- 
teraction if its eccentricity is not small to begin with. This would 
be the case for a planet at the end of the chaotic phase. 

Additional simulations are undergoing in order to perform a 
deeper exploration of the parameter space to identify the initial 
configurations of the planets and disk that lead to planet-planet 
scattering rather than to a stable resonant trapping. 

Acknowledgements. We thank Phil Armitage for his useful comments and sug- 
gestions as referee of the paper. 
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Fig. 1. Orbital evolution of the three planets while embedded in the disk. The top plots show a model where Zq = 0.4Emmsn and 
fli = 3.0AU, 02 - 6.7AU, (23 = 10.5AU The middle plots show a case with Zq = 0.65Emmsn and oi - 2.0AU, 02 = 6.5AU, as = 
lO.OAU The bottom plots illustrate the evolution when Xq - Smmsn and oi = 3.0AU, 02 - 5.5AU, 03 - lO.OAU. 



